Stability of spherically symmetric solutions in modified theories of gravity 
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In recent years, a number of alternative theories of gravity have been proposed as possible res- 
olutions of certain cosmological problems or as toy models for possible but heretofore unobserved 
effects. However, the implications of such theories for the stability structures such as stars have 
not been fully investigated. We use our ''generalized variational principle'' , described in a previous 
work pj, to analyze the stability of static spherically symmetric solutions to spherically symmetric 
perturbations in three such alternative theories: Carroll et aJ.'s f(R) gravity, Jacobson & Mat- 
tingly's "Einstein-aether theory", and Bekenstein's TeVeS. We find that in the presence of matter, 
f(R) gravity is highly unstable; that the stability conditions for spherically symmetric curved vac- 
uum Einstein-aether backgrounds are the same as those for linearized stability about flat spacetime, 
with one exceptional case; and that the "kinetic terms" of vacuum TeVeS are indefinite in a curved 
background, leading to an instability. 

PACS numbers: 04.20.Fy, 04.40.Dg, 04.50.+h 



I. INTRODUCTION 

The idea of using a variational principle to bound the 
spectrum of an operator is familiar to anyone who has 
taken an undergraduate course in quantum mechanics: 
given a Hamiltonian H whose spectrum is bounded be- 
low by Eq, we must have for any normalized state \ip) 
in our Hilbert space Eq < (ijj\H\ip). We can thus obtain 
an upper bound on the ground state energy of the sys- 
tem by plugging in various "test functions" \tp), allowing 
these to vary, and finding out how low we can make the 
expectation value of H. 

This technique is not particular to quantum mechan- 
ics; rather, it is a statement about the properties of self- 
adjoint operators on a Hilbert space. In particular, we 
can make a similar statement in the context of linear field 
theories. Suppose we have a second-order linear field the- 
ory, dependent on some background fields and some dy- 
namical fields ip a , whose equations of motion can be put 
into the form 



dt 2 



(1) 



for some linear spatial differential operator T. Suppose, 
further, that we can find an inner product (•, •) such that 
T is self-adjoint under this inner product, i.e., 



(2) 



for all ipi and ip2- Then the spectrum of T will corre- 
spond to the squares of the frequencies of oscillation of 
the system; and we can obtain information about the 
lower bound of this spectrum via a variational principle 
of the form 



, )2< MM) 



(3) 
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Moreover, if we find that we can make this quantity nega- 
tive for some test function ijj, then the system is unstable 
0; and we can obtain an upper bound on the timescale 
of this instability by finding how negative Uq must be. 

One might think that this would be a simple way to 
analyze any linear (or linearized) field theory; in prac- 
tice, however, it is easier said than done. A serious prob- 
lem arises when we consider linearizations of covariant 
field theories. In such theories, the linearized equations 
of motion are not all of the form ([1]); instead, we ob- 
tain non-deterministic equations (due to gauge freedom) 
and equations which are not evolution equations but in- 
stead relate the initial data for the fields tp a (since co- 
variant field theories on a spacetime become constrained 
when decomposed into "space + time".) Further, even 
if we have equations of the form it is not always 
evident how to find an inner product under which the 
time-evolution operator T is self-adjoint (or even if such 
an inner product exists.) Thus, it would seem that this 
method is of limited utility in the context of linearized 
field theories. 

In spite of these difficulties, Chandrasekhar success- 
fully derived a variational principle to analyse the stabil- 
ity of spherically symmetric solutions of Einstein gravity 
with perfect-fluid sources [H, 0|- The methods used in 
these works to eliminate the constraint equations, put 
the equations in the form ([T]), and obtain an inner prod- 
uct seemed rather particular to the theory he was ex- 
amining, and it was far from certain that they would 
eneralize to an arbitrary field theory. In a recent paper 
l|, we showed that these methods did in a certain sense 
"have to work out" , by describing a straightforward pro- 
cedure by which the gauge could be fixed, the constraints 
could be solved, and an inner product could be obtained 
under which the resulting time-evolution operator T is 
sclf-adjoint. We review this procedure in Section|TT]of the 
paper. We then use this "generalized variational princi- 
ple" to analyze three alternative theories of gravity which 
have garnered some attention in recent years: f(R) grav- 
ity, the current interest in which was primarily inspired 
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by the work of Carroll et al. [|[; Einstein-aether theory, 
a toy model of Lorentz-symmetry breaking proposed by 
Jacobson and Mattingly 0, 0|; and TeVeS, a covariant 
theory of MOND proposed by Bekenstein [8| . We apply 
our techniques to these theories in Sections IIII1 IIV1 and 
IVl respectively. 

We will use the sign conventions of [j| throughout. 
Units will be those in which c = G = 1. 



(i.e., ebcd — n a e a bcd), this can be written in terms of to a 
instead: 

fi=-f(w\)e. (9) 

In performing the calculations which follow, it is this sec- 
ond expression for f2 which will be most useful to us. 



II. REVIEW OF THE GENERALIZED 
VARIATIONAL PRINCIPLE 

A. Symplectic dynamics 

We first introduce some necessary concepts and nota- 
tion. Consider a covariant field theory with an action of 
the form 




where is a scalar depending on some set 'J of dynam- 
ical tensor fields including the spacetime metric. (For 
convenience, we will describe the gravitational degrees of 
freedom using the inverse metric g ah rather than the met- 
ric g a b itself.) To obtain the equations of motion for the 
dynamical fields, we take the variation of the four- form 
Ce with respect to the dynamical fields ^: 

5{Ce) = + V a 6» Q [*,5*])e (5) 

where a sum over all fields comprising ^> is implicit in 
the first term. Requiring that SS = under this varia- 
tion then implies that the quantities £^ vanish for each 
dynamical field '5. 

The second term in (|S|) defines the vector field 
6 a [^f, 5^]. The three-form 8 dual to this vector field 
(i.e., Obcd = 6 a £abcd) is the "symplectic potential cur- 
rent." Taking the antisymmetrized second variation of 
this quantity, we then obtain the symplectic current 
three-form u) for the theory: 

Si^f, tf 2 *] = Si6[^, <S 2 *] - <5 2 0[#, <*i#]- (6) 

In terms of the vector field u> a dual to u>, this can also be 
written as 

W a £ a bcd = Sl{02 e abcd) ~ ^(Ox^abcd)- (7) 

The symplectic form for the theory is then obtained by 
integrating the pullback of this three-form over a space- 
like three-surface E: 

n[\Mi\M 2 *] = / w[*;tfi#,<f2*]- (8) 

If we define n a as the future-directed timelike normal 
to E and e to be the induced volume three-form on E 



B. Obtaining a variational principle 

In we presented a procedure by which a varia- 
tional principle for spherically symmetric perturbations 
of static, spherically symmetric spacetimes could gener- 
ally be obtained. For our purposes, we will outline this 
method; further details can be found in the original pa- 
per. 

The method described in [l| consists of the following 
steps: 

1. Vary the action to obtain the equations of motion 
(£a)ab corresponding to the variation of the met- 
ric, as well as any other equations of motion £a 
corresponding to the variations of any matter fields 
present. This variation will also yield the dual 6 a 
of the symplectic current potential; take the anti- 
symmetrized variation of this quantity (as in ([7])) 
to obtain the symplectic form ([9]). 

2. Fix the gauge for the metric, and choose an ap- 
propriate set of spacetime functions to describe 
the matter fields. Throughout this paper, we will 
choose our coordinates such that the metric takes 
the form 

ds 2 = _ e 2*(r,t) dt 2 + e 2A(r,t) ^2 + ^2 (1Q) 

for some functions $ and A. 1 Our spacetimes will 
be static at zero-order (i.e., {d/dt) a is a Killing vec- 
tor field at zero-order), but non-static in first-order 
perturbations. We will therefore have $(r, t) — 
$(r) + 0(r, t), where <f> is a first-order quantity; sim- 
ilarly, we define the first-order quantity A such that 
A(r,i) = A(r) + A(r, t). In other words, 

8g tt = 2a- 2 *W<j>(r,t) (11) 

and 

Sg rr = -2e- 2A ^X(r,t). (12) 

Similarly, all matter fields will be static in the back- 
ground, but possibly time-dependent at first order. 



Such a set of coordinates can always be found for a spherically 
symmetric spacetime [Toll , 
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3. Write the linearized equations of motion and the 
symplectic form in terms of these perturbational 
fields (metric and matter.) 

4. Solve the linearized constraints. One of the main 
results of [1] was to show that this can be done 
quite generally for spherically symmetric perturba- 
tions off of a static, spherically symmetric back- 
ground. Specifically, suppose the field content 

of our theory consists of the inverse metric g ab and 
a single tensor field A ai - " an b 1 ...b m - (The general- 
ization to multiple tensor fields is straightforward.) 
Let (£U)ai...a„ " bm denote the equation of motion 
associated with A ax '" a " n bx...b m - We define the con- 
straint tensor as 

C c d — 1 2{£g)ccI 

— 9ce 22 '"" // ..../.../. . {£A)a 1 ...a n bl "' e "' hm 

i 

■ H,.y,.V' : - ' " /, ....... (£ A )a (13) 

i 

where the summations run over all possible index 
"slots", from 1 to n and from 1 to m for the first 
and second summation respectively. It can then be 
shown that if the background equations of motion 
hold, and the matter equation of motion also holds 
to first order, the perturbations of the tensor C a b 
will satisfy 

<>L = -r 2 e*- A 5C rt (14) 

and 

^ = -r 2 e A ^5C tt (15) 
or 

for some quantity F which is linear in the first- 
order fields. Moreover, the first-order constraint 
equations SCtt — 5Ct r — will be satisfied if and 
only if F = 0. We can then solve this equation 
algebraically for one of our perturbational fields, 
usually the metric perturbation A. We will refer to 
the equation F — as the preconstraint equation. 

5. Eliminate the metric perturbation <fi from the equa- 
tions. As <f> cannot appear without a radial deriva- 
tive (due to residual gauge freedom) 2 , we must find 
an algebraic equation for d<j)/dr. The first-order 
equation 6C rr = 0, will, in general, serve this pur- 
pose |l[. Use the above relations for A and dcj)/dr 



to eliminate the metric perturbations completely 
from the perturbational equations of motion and 
the symplectic form, leaving a "reduced" set of 
equations of motion and a "reduced" symplectic 
form solely in terms of the matter variables. 

6. Determine if the reduced equations of motion take 
the form |T]). If so, read off the time-evolution op- 
erator T. 

7. Determine if the symplectic form, written in terms 
of the reduced dynamical variables ip a , is of the 
form 

0(*;^,^) = ^W^(^-^) (16) 

for some three- form W Q| g. Then (as we showed in 
[l|) we must have W ^ = and 

/ W Q/ #«TV 2 7 = / W^T'V/. (17) 

JE JE 

8. Determine whether W a ^ is positive definite, in the 
sense that 

/ W a/ ^ V > (18) 
Je 

for all ip a , with equality holding only when ip a = 
0. 3 If this is the case, we can define an inner prod- 
uct on the space of all reduced fields: 

{ih,ih) = J Wc^t^. (i9) 

Equation (fT7|) then shows that T is a symmetric 
operator under this inner product. Thus, we can 
write down our variational principle of the form 
©■ 

It is important to note that each step of this procedure 
is clearly delineated. While the procedure can fail at cer- 
tain steps (the reduced equations of motion can fail to 
be of the form ([T]), for example, or W Q ^ can fail to be 
positive definite), there is not any "art" required to apply 
this procedure to an arbitrary covariant field theory. (In 
practice, as we shall see, there are certain shortcuts that 
may arise which we can exploit; however, the "long way" 
we have described here will still work.) In the follow- 
ing three sections, we will use this formalism to analyze 
the stability of f(R) gravity, Einstein-aether theory, and 
TeVeS. 



2 The exception to this statement is when the matter fields under 
consideration have non-zero components in the i-direction. In 
this case, <j> and its derivatives can appear in linear combination 
with perturbations of the matter fields; however, a redefinition 
of the matter fields will suffice to eliminate such cases. 



If W a( 3 is negative definite in this sense, we use the negative of 
this quantity as our inner product, and the construction proceeds 
identically. 
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III. f(R) GRAVITY 
A. Theory 

In f(R) gravity, the Ricci scalar R in the Einstein- 
Hilbert action is replaced by an arbitrary function of R, 
leaving the rest of the action unchanged; in other words, 
the Lagrangian four-form C is of the form 



1 

16^ 



f(R)e + C mat [A,g ab ] 



(20) 



where A denotes the collection of matter fields, with ten- 
sor indices suppressed. Taking the variation of this ac- 
tion with respect to the metric, we obtain the equation 
of motion 

f'(R)R ab - l -f{R)g ab 

- V a V 6 /'(fl) + gabOf'(R) = 8TrT ab (21) 

where T ofc , given by 



<5£mat = ~\{T ab 5g ab )e, 



(22) 



is the matter stress-energy tensor, and f'(R) = df/dR. 

This equation is fourth-order in the metric, and as such 
is somewhat difficult to deal with. We can reduce this 
fourth-order equation to two second-order equations us- 
ing an equivalent scalar-tensor theory [ll], [l2|]. 4 This 
equivalent theory contains two dynamical gravitational 
variables, the inverse metric g ab and a scalar field a, in 
addition to the matter fields. The Lagrangian is given by 



£ = tt - {f'i a )R + f( a ) - a f'( a )) e - 



C mat [A,g ab ]. 



(23) 

Varying the gravitational part of action with respect to 
g ab and a gives us 

SC = {(£ G )abSg ab + £ a Sa + E A 5A + V a 9 a ) e (24) 

where £a denotes the matter equations of motion, 



1 

167T 



f'(a)G ab - V„V 6 /'(a) + 9abOf'(a) 
1 



75ab(/(") - "/'(")) - 87rT afc 



(25) 



and the vector 9 a is our symplectic potential current: 

0° = /'(«)^Ei„ + Cat 

+ ^ {(Vbf(a))5g ab - (\7 a f'(a))g bc Sg bc ) . (27) 

The vector #^ lat above is the symplectic potential current 
resulting from variation of £ ma t, and 9^ in is the symplec- 
tic potential current for pure Einstein gravity, i.e., 



7 Ein 



1 

lOTT 



{ 9 bcV a Sg bc - V b 5g ab ) 



(28) 



The equations of motion are then given by (£a)a b = 
and £ a = 0. Assuming that /"(a) ^ 0, this second equa- 
tion implies that R = a, and substituting this relation 
into (pj)) yields the equation of motion obtained in ([2"Tj) . 
Hereafter, we will use the form of the equations obtained 
from the action (l23l) . 



and 



B. Obtaining a variational principle 

Before we examine the stability of spherically symmet- 
ric static solutions in f(R) gravity with perfect fluid mat- 
ter, we must first consider the question of whether phys- 
ically realistic solutions exist. In particular, do there ex- 
ist solutions to the equations of motion which reproduce 
Newtonian gravity, up to small relativistic corrections? 
There has recently been a good deal of debate on this 
subject [3 [H, [Ti, [T3, HI. We will not comment di- 
rectly on this controversy here except to say that if f(R) 
gravity (or any theory) does not allow interior solutions 
with R w 8nGp to be matched to exterior solutions with 
R close to zero, then it is difficult to see how such a the- 
ory could reproduce Newtonian dynamics in a nearly-flat 
spacetime. We will therefore give the theory the "benefit 
of the doubt," and assume that such solutions exist. 

An initial attempt to address the stability of such solu- 
tions was made by Dolgov and Kawasaki [19(; their per- 
turbation analysis implied that stars would be extremely 
unstable in Carroll et ai.'s f(R) gravity, with a charac- 
teristic time scale of approximately 10 -26 seconds. Their 
results, while suggestive, nevertheless failed to take into 
account the constrained nature of the theory: the stress- 
energy tensor, the metric, and the scalar a cannot all 
be varied independently. The imposition of a constraint 
can, of course, change whether or not a given system is 
stable; as a trivial example, consider a particle moving 
in the potential V(x, y) = x 2 — y 2 , with and without the 
constraint y = const. 5 In what follows, we will show 
that Dolgov and Kawasaki's conclusion is, nonetheless, 



£ a = f"(a)(R-a), 



(26) 



4 This equivalence was also used in the special case f(R) = R 
2A + aR 2 in ffj]. 



5 In certain scalar-tensor theories |20H . the scalar perturbations 
decouple from the metric and matter perturbations; it is then 
legitimate in such theories to "ignore" the constraints if we con- 
cern ourselves only with the scalar field. However, this de- 
coupling does not occur in f(R) gravity theory. To see this, 
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correct: stars in CDTT f(R) gravity do in fact have an 
ultra-short timescale instability. 

To describe the fluid matter, we will use the "La- 
grangian coordinate" formalism, as in Section V of 
In this formalism, the fluid is described by considering 
the manifold M. of all fluid worldlines in the spacetime, 
equipped with a volume three-form N. If we introduce 
three "fluid coordinates" X A on Ai, with A running 
from one to three, then the motion of the fluid in our 
spacetime manifold M is completely described by a map 
X ■ M — > M. associating with every spacetime event x 
the fluid worldline X A (x) passing through it. The mat- 
ter Lagrangian is then given by 

C = -g{v)e (30) 

where v, the "number density" of the fluid, is given by 



<*«'■■■ if + 



2<9$ 1 
r dr r 



Z^~^ 2A -!))/'(«) 



+ -e 2A (a/'(a) - f(a)) - e 2A {g'u - g) = 0, (34b) 



G) 



2 -2A 

r e 



f-f(a) + (?* - ^ + I) JLf'fa) 
dr 2 \ dr dr r ) dr 



<9 2 $ <9$ OA 
dr 2 dr dr 

*.2 



<9$ 

dr 



1 [d<5> <9A N 



/'(«) 



r \ dr dr 
(af'(a) - /(a)) - r 2 {g'v -g) = 0, (34c) 



v 2 = \N abc N a 
6 



(31) 



In turn, N abc , the "number current" of the fluid, is given 
by 



N abc = AT ABC (^)V X /1 V 6 X a V c X . (32) 

We will be purely concerned with spherically symmetric 
solutions and radial perturbations of the fluid; thus, we 
will take our Lagrangian coordinates to be of the form 

X R = r, X B = (9, X* = ip (33) 

in the background, and consider only perturbations 
SX R ee £(r,t) at first order (i.e., SX e = SX* = 0.) 

We can then easily obtain the background equations of 
motion; these are 



G)tt 



2$-2A 



2dA 1 

r dr r 



2 (e 2A -!))/» 



ie 2 *(a/'(a)-/(a))-e 2 *^0, (34a) 



we can rewrite the Lagrangian in terms of a rescaled metric 
gab = n- 2 (cr) 9 a6 = e a g al > , putting it m tnc form [H],[Ii 

C=( -±-R - iv a <rV a CT - V(a)) e+C aia , t [A,e' r / 2 ^g ab ] (29) 



V16tt 



where a is related to the scalar a by /'(a) = e"^ 2 ^^ and the 
exact form of the potential Via) is determined by the function 
f(R)- The decoupling of the scalar perturbations, however, re- 
quires V(ao) = and Q'(ao)/Q = 0, where co is the background 
value of a. The former condition will not hold for a generic 
f{R), and the latter condition will not hold for any f(R) for 
which f"(R) ^ 0. 



£ a = 2/»e 



-2A 



9 2 $ d<S> dA 

dr 2 dr dr 



- f"(a)a = 0, 



and 



(fx)fl = g' 



.du d* , 

-5- + -z-Q v = a, 

dr dr 



(34d) 



(34e) 



where g' = dg/dv. These equations are not all inde- 
pendent; in particular, the Bianchi identity implies that 
Q34cj) is automatically satisfied if the other four equations 
are satisfied. Note that under the substitutions g — > p, 
g'v — g — > P, the "matter terms" in these equations take 
on their familiar forms for a perfect fluid. 

We next obtain the symplectic form for the theory. 
The contribution #g rav to the symplectic potential cur- 
rent from the gravitational portion of the action is given 
by ([23; to obtain the symplectic form, we then take the 
antisymmctrized second variation of 6, as in l[7]). Per- 
forming this variation, we find that 

+ [5 1 (/'(a))V^ 2ff cd - V fc (*i(/'(a))) S 2 g cd ] 

x (g ab g cd - 6 a c S b d ) - [1 <-» 2] (35) 

where Wg in is defined, analogously to 0g in , to be the 
symplectic current associated with pure Einstein grav- 
ity. This is equal to [H| 

< n = ^ a 6c d e /(^5 &c V^ig e/ - 8 ig bc V d 5 2 g e f ), (36) 

where 



he ef = yS a e S c g bf - -g a g be g cf 

■-5 a b S d c g e f - ^S a eS d f g bc + ^g ad gb C gen ■ (37) 
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The contribution to the symplectic current from the 
matter terms in the Lagrangian coordinate formalism 
was calculated in [l[; we simply cite the result for the 
t-component of oj a in a static background here: 

i Q < at = -t a ^-N ABC V b X B V c X c [S l9 ad S 2 X A N d bc 

+35 2 X A V^ a (N DEF 6 1 X D )V b X E V^X F - [1 <-► 2] 

"(38) 

where the antisymmetrization in the second term is over 
the tensor indices only (not the fluid-space indices). 
Writing out u> in terms of our perturbational variables, 
we have 



56 

at 



6 



96 
dt 



6 



-2# 



<9A 2 dbi 
~di dt 



dXi db 



+ ^Xi (39) 



dt dt 



where we have defined b = S(f'(a)) = f"(a)6a. (Note 
that the t-component of the symplectic current about 
a spherically symmetric static solution vanishes in pure 



Einstein gravity, i.e., w 



Ein 



0. 



The first main step towards obtaining a variational 
principle is to solve the linearized constraints. To do this, 
we must calculate the constraint tensor C ab , as defined 
in (flU|) . Since all the fields in our theory other than the 
metric are scalars with respect to the spacetime metric, 
the resulting expression is particularly simple: 

C ab = 2(£ G )afc, (40) 

where (£c)ab, the metric equation of motion, is given by 
([25]). We can then find algebraic equations for A and 
dcf>/dr. The algebraic equation for A will be given by the 
solution of the equation F = 0, where F is given by (jT5J) 
and (|15[) . Similarly, an algebraic equation for d(j)/dr can 
be found by solving the equation 

SC rr = (S£ G )„ = 0. (41) 

The situation is thus very similar to the example of 
pure Einstein gravity coupled to perfect fluid matter de- 
scribed in [l| ; only the precise form of the perturbational 
equations of motion F — and (5£a) rr — are different. 
In the f(R) gravity case, these equations are 

9h M-b-f^M (42) 



F = S\ 



dr dr 



where we have introduced the new quantity 



(44) 



for notational convenience. Equation (|42f then implies 
that our algebraic equation for A is 



x = s _ x (db_f b + 8ire .A M 
\ or or 



(45) 



We could plug this result into P3"| to obtain an alge- 
braic equation for d<p/dr\ however, the resulting expres- 
sion is somewhat complicated. In fact, we can derive a 
simpler expression for d<j)/dr. To do so, we note that 
(£o)ab - \g a b(£G)c C = if and only if {£ G )ab = 0. In the 
case of f{R) gravity, this equation is given by 



f'(a)R ab - -g ab nf(a) - V.V 6 /'(a) 

+ \gab (f(a) - f'(a)a) - Sir (r ab - X -g ah T c 



= 0. 

(46) 



Using the trace of this equation to eliminate the □/'(a) 
term, along with the equation R = a, we can show that 



f'(a)R ab + -g ab (/(a) - 2/'(a)a) - V„V 6 /'(a) 
6 



I* ( T a b - ^g ab T c c ) = 0. (47) 



To zero-order, the ##-component of this equation is 



and 



(S£ G )r 



-2A 



-2A 



S 7T- e 
or 



5$ 



-2<1> 



-2A 



2 9$ 



(Pb 
dt 2 



-2A 



d<f> 2\ db 



dr r I dr 



i 



-2A\ 



",.2 



r dr r 2 
d£ fdk 



di 



di 



2 1 dv\ 
r v dr ) 



2 



(43) 



V(a)(^-^ + i(e- 
r \ dr or r 



1 



1 d 



G 



and to first order, it is 



r dr" 



87T 9A 



(48) 
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1 , , fox d$ , 2 2A 

r \ or or r 



dA 
dr 



<9$ 

9r 



(e 2A - 1) - ie 2A 



/"(a) 



2a 



, 1 36 



1 



3 6 



2A 



(/(«) 



Tafia) ~ 16.,) A - y ^ f^ + (_ + - + -_ U -A 



0. (49) 



This equation, with A given by (|4"5"|) . can then be solved 
algebraically for d(j)/dr in terms of the matter variables 
and their derivatives. 

Our next step is to obtain the reduced matter equa- 
tions of motion, i.e., eliminate the metric degrees of free- 
dom A and d<j)/dr from the remaining first-order equa- 
tions of motion. These remaining equations of motion 
are the "scalar" equation R — a = 0, which at first order 



d 2 4> d<S>d\ dAdcj) 2 <9$ dt 
dr 2 dr dr dr dr dr dr ' r V dr dr 



, 2 2A1 2A-2*^ 2 ^ 



e 2A | aX + 



2 ^d,\ On 

b 



2f"{a) 



= 



(50) 



and the matter equation of motion, which as in the case 
of pure Einstein gravity is 



„2A-2$ d 2 £ 

dt 2 



d <9$ 

dr dr 



Q v 



dr 

d A 

dr 



+ 



1 dv 
v dr 



dA 
dr 



0. 



A 
(51) 



We can then eliminate the gravitational equations of mo- 
tion using (|45|) and l[49|) . leaving equations solely in terms 
of the matter variables b and £. Equivalcntly, we can 
write our "reduced" equations in terms of b and a new 
variable £, defined by 



(52) 



(The utility of this new variable will become evident 
when we obtain the reduced symplectic form.) Perform- 
ing these algebraic manipulations, the resulting reduced 
evolution equation for £ takes the form 



At 



d^C 

dr 2 



dr 



db 
dr 



- A 2 ^r L + AzC + A 4 ^ + A 5 b 

(53) 

where the coefficients Ai are dependent on the back- 
ground fields. Similarly, the evolution equation for b 
takes the form 



-2* 



d% 



Bi 



d 2 b 
dr 2 



B 



db 
dr 



B 3 b + a 



d( 
dr 



B 5 (. (54) 



We see from the form of the above equations that f(R) 
gravity has cleared another hurdle required to have a 



valid variational principle: the reduced equations do in 
fact take the form (TTJ) , containing only second derivatives 
in time and only up to second radial derivatives. 

We can also apply these constraint equations to the 
symplectic form. For the symplectic form, we only need 
to eliminate A using the first-order constraint equation 
(I45p . Performing an integration by parts to eliminate 
the mixed derivatives that result, and applying the back- 
ground equations of motion, we obtain our reduced sym- 
plectic form: 



= 4tt / drr 2 e A+ * 



,-2* c-2 



dt 



dt 

-[!<-> 2]. (55) 



As noted above (|16p. this reduced symplectic form de- 
fines a three- form W a /3- For a valid variational principle 
to exist, this W Q ^ must be positive definite in the sense 
of (fl"8|) . We can see that the W Q/ 3 defined by (|55|) is 
positive definite if and only if f'(a) > and g'v > 
in our background solutions. The latter condition will 
hold for any matter satisfying the null energy condition, 
since q'v = p + P; however, the former condition must 
be checked in order to determine whether a valid varia- 
tional principle exists. 6 For the particular f(R) chosen 
by Carroll et al. [|, we have f{R) = 1 + ///i? 2 > 0, and 
so W Q £j is always positive definite in the required sense 
for this choice of f(R). 

All that remains is to actually write down the varia- 
tional principle for f(R) gravity. As noted above ([3]), the 
variational principle will take the form 



2 . OW) 
Wn < 



(56) 



where T is the time-evolution operator. In our case, ip 
denotes a two-component vector, whose components are 
the functions £ and b. The inner product in which T 
is self-adjoint can be read off from {55]), and thus the 



6 Should /'(«) fail to be positive in the background solutions, all is 
not necessarily lost; we can still attempt to analyze the reduced 
equations that we have obtained. See Section [V] for an example 
of such an analysis. 
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denominator of ([56)1 will be 

O,V0 = 4tt / drr 2 e A -' s ' 



we find that 



(57) 

The numerator of (f56|) is, of course, rather more compli- 
cated. After multiple integrations by parts and applica- 
tions of the background equations of motion, we can put 
this quantity in the form 



(tp, Ttj}) = 47r J dr 
8b 



or 



+e 3 c 



4^) +C 4 6 2 +C 5 C 2 + 2(',k 



(58) 



where the Ci coefficients are functions of the background 
fields. The exact forms of the Ci are given in Appendix 

m 



C. Discussion 

All of our results thus far have been independent of 
the choice of f(R) (assuming, of course, that f"(R) ^ 0.) 
In the case of Carroll et al.'s f(R) gravity, we can now 
use this variational principle to show that this theory is 
highly unstable for Newtonian solutions. For the choice 
f(R) = R - p 4 /R, we have f(R) = 1 + p 4 /R 2 and 
f"(R) = — 2p 4 /R 3 . Moreover, for a quasi-Newtonian 
stellar interior, we will have a = R tts p, where p is the 
matter density in the star. In particular, this implies that 



C 4 w 2e* +A 5 



? - 2 (f(«)) 2 „ «+Ao-a^ fw1 



since the above term will dominate over all the others in 
C4. (Note that for a star with the density of the Sun, 
and the choice of p « 10 -27 m _1 made by Carroll et al, 
p/p 2 ~ 10 10 .) This implies that for a test function ip 
with £ set to zero, we will have 



< 



Jdr 



„A+* 



s- 



P 7) 2 



6e 



-2A 



6 / dr e A -^S- 2 b 2 



(60) 



(note that f'{a) — 1 + p i /p 2 ~ 1 in the stellar inte- 
rior.) As a representative mass distribution, we take the 
Newtonian mass profile of an n = 1 polytrope: 



P{r) = pa 



(61) 



where i? is the radius of the star and po is its central 
density. We take p to be of a typical stellar density, 
Po ~ 10 -24 metres -2 . Numerically integrating (|59p with 
a test function of the form 



1-fc r<i? 
r > i? ' 



(62) 



w 2 < -9 x 10 35 metres" 



(63) 



which corresponds to an instability timescale of t k 4 x 
10 -27 seconds. This timescale is of the same magnitude 
as the instability found in fl9j | . 

We see, then, that for the f(R) originally chosen by 
Carroll et al., the theory is extremely unstable in the 
presence of matter. Moreover, a similar argument will 
obtain for any choice of f(R) for which quasi-Newtonian 
solutions exist and for which f"(R) is sufficiently small 
and negative at stellar-density curvature scales. We can 
always pick a test function b(r) lying entirely inside the 
stellar interior such that db/dr is of order b/R. Thus, 
if \f'(p)/f"(p)\ > 1/i? 2 for a typical stellar density p, 
the b 2 term in (fS"9"|) will dominate the (db/dr) 2 term, and 
the resulting lower bound on w 2 will then be of order 
/'(a) //"(a). If the choice of /(a) results in this quantity 
being negative, quasi-Newtonian stellar solutions will be 
unstable in the corresponding theory. We can thus rule 
out any theory (e.g., [aHH) for which f"(p) is sufficiently 
small and negative. 

This result seems to be closely related, if not identical 
to, the "high-curvature" instabilities found in cosmologi- 
cal solutions by Song, Hu, and Sawicki [23j. Particularly 
telling is the fact that the instability time scale found in 
their work is proportional to y/\f"(R)/f'(R)\, the same 
time-scale found in the present work. 



IV. EINSTEIN-jETHER theory 



A. Theory 

Einstein-aether theory 0, 0] was first formulated as a 
toy model of a gravitational theory in which Lorentz sym- 
metry is dynamically broken. This theory contains, along 
with the metric g a b, a vector field u a which is constrained 
(via a Lagrange multiplier Q) to be unit and timelike. 
The Lagrangian four-form for this theory is 

C = (~^R + K ab cd V a u c V b u d + Q(u a u a + 1)^ e 

+ C mat [A,g ab ,u a ] (64) 

where A denotes any matter fields present in the theory, 
and 

K ab cd = c 1 g ab g cd +c 2 S a c S b d + c 3 S a d S b c -c i u a u b g cd . (65) 

The Ci constants determine the strength of the vector 
field's coupling to gravity, as well as its dynamics. 7 (Note 



7 Note that our definitions of the coefficients Ci differ from those in 
[tj Fil [24l. |25(| . as do the respective metric signature conventions. 
The net result is that to translate between our results and those 
of the above papers, one must flip the signs of all four coefficients. 
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that in the case C3 = —c\ > and C2 = C4 = 0, we have 
the conventional kinetic term for a Maxwell field; this 
special case was examined in 26] , prior to Jacobson and 
Mattingly's work.) In the present work, we will work in 
the "vacuum theory", i.e., in the absence of matter fields 
A. 

Performing the variation of the Lagrangian four-form, 
we find that 

6C = ((£ G )a b Sg ab + (£ u ) a 6u a + £ Q SQ + V a 6 a ) e (66) 
where 

(£c)ab = J^Gab + V c (j C ( a Ub) + J(ab)U C - 

+ Ci (V a M C V b M c - V C U a V c U b ) + CiilaUb 

- Qu a u b - ^gabJ c d^cU d , (67a) 



(£ u ) a = -2V b J b a - 2c 4 u b V a u b + 2Qu a , (67b) 



(£ Q ) = u a u a + 1, 



(67c) 



and 



9° in +2J a b 8u b +(J b a u c - J\u c - J bc u a ) 5g bc . (68) 



In the above, we have introduced the notation u a = 
u b V b u a and J a b = K ac bd V c u d ; 9% in is defined by ([28]). 
as above. Equation (|67cp is the constraint that u a is unit 
and timelike, while (|67a[) and (|67b[) are the equations 
of motion for g ab and u a , respectively. If desired, we can 
eliminate the Lagrange multiplier Q from these equations 
by contracting (|67b[) with u a , resulting in the equation 



= -u a V fc J\ - c 4 u a ii b V a u b 



(69) 



We now take the variation of the symplectic potential 
current to obtain the symplectic current. The resulting 
expression can be written in the form 



^Ein 



(70) 



where 



y Ein 



is the usual symplectic current for pure Ein- 



stein gravity (given by (p?6"|) ), and wj cc is obtained by tak- 
ing the antisymmetrized variation of the last two terms 
in ([68]) . The exact form of this tensor expression is rather 
complicated, and can be found in Appendix [Bj 



i.e., the first-order perturbation to u* is — </>e~* while 
the first-order perturbation to u r is v. (Note that un- 
der these conventions, g ab u a u b = — 1 to first order, as 
required.) The background equations of motion can then 
be calculated to be 




-1) 
dr V dr 



and 



(£a)r 



2 9$ 
r dr 



1 



-,2A 



1 



C14 (d® 
~2~ Vdr 



(72) 



(73) 



where we have defined C14 = C1+C4 and used the equation 
5$ / 8A 2" 



,-2A 



C3 



dr \ d 



-(ci + c 3 + 2c 4 ) ( — J -c 3 ^j 



(74) 



to eliminate Q. The final equation of motion, (£ u ) a = 0, 
is satisfied trivially if the aether is static. 

In terms of these variables, the i-component of ui a can 
be calculated to be 



uj 1 = 2e 



-2<I> 



Cl23^T- A 2 + Ci23e 



d\i 
~dT 



dr 



dA . .8$ 2\ 

e~ I Cl23^- + (Cl23 - c 14)"5- + C 2 - Vi\ 2 

dr dr r 



dv 2 

+ci 4 e vi— he c u vi-— 

dr dt 



„2A, 



[1^2] (75) 



where we have defined C123 = C1+C2+C3. In what follows, 
we will assume that C14 ^ 0, and, except where otherwise 
noted, that C123 7^ as well. 

Our next step, as usual, is to solve the constraints. 
However, in this theory we have the added complication 
of the presence of a vector field. This means, in particu- 
lar, that the tensor C ab (as defined in ([13]) ) is not merely 
proportional to (£c)ab, as in the previous section, but is 
instead 



B. Obtaining a variational principle 

As before, we are primarily concerned with the t- 
component of uj a , and its form in terms of the pertur- 
bational variables. We will use the usual spherical gauge 
(flOf for g ab . We will further assume that u a oc t a in 
the background solution (the so-called "static aether" as- 
sumption), and that to first order we have 

u" = (e~* -<£e-*,u,0,0), (71) 



C ab = 2(£ G )ab + U a (£u)b (76) 

= ^ Gab + 2Vc ( J (ab) uC ~ J(a C U b) ) 
+ 2 (V c J C [ a U b ] + V C M( Q J C b) ) 

+ 2ci {V a u c V b u c - V c u a V c u b ) 

- 2c 4 (u a ii c V b u c - u a u b ) - g a bJ c <N c u d . (77) 

Writing out SC r t in terms of our perturbational quan- 
tities, we find that it is indeed a total time derivative, 
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with {SEcjrr = and (S£ u ) r = 0; in terms of the pertur- 

bational variables, these are 

or ) o pfi\ 

8l) 8rh\ [OOGjrr = ,e A + Ci 23 e —-r 

C14-5- 7T + c i23e 



r or J or drdt 



Note that this quantity depends on 0; and its deriva- 
tives can appear in the preconstraint equation F = _|_ e 2A-* | ^£2. Cl23 £^. _(_ ( Cl23 _ C14) ^ J — (79) 



when our background tensor fields have non-vanishing \ r dr dr J dt 

i-componcnts, as noted in Footnote^ 

The remaining non-trivial equations of motion are and 



-{8£ u ) r = c u e 



2A-2* 



dt 2 



C123 



d 2 v 



dr 2 

2 <9A 

C123 I - + -5- 
r or 



Cl23e 



d 2 \ 



drdt 
<9$\ dv 
dr J dr 



■ cue 



drdt 



(C 1 23-C 14 )--( C1 +C 3 )-1- 



x 9 2 $ d 2 A dAd<i> if. dA ,<9$ 

(C123 " C 14)^- + C 123^J + c 14 — — + I ( ( C l + C 3)— + ( C 3 ~ Ci) — 



dr dr 
V 



dr 



dr 



2 

Cl23 ^2 



v. (80) 



While we could follow the methods outlined in Section 
[XT] to reduce these equations to the basic form |T]), it is 
actually simpler to pursue a different path. If wc solve 
([75)1 for d(j>/8r, rather than A as usual, and plug the 
resulting expressions into (l79|) and (j80|) . there result the 
equations 



dt 



„2$-2A 



2 
Cl4 



-tt — + 1 A = 



(81) 



and 



dip 
dr 



ci + c 3 + 
C123 
(ci- 



1 2 

r 

" C3 



dr 

1)(C123 



2c 2 



C123 



-v = 0, (82) 



where we have introduced the new variable ip, defined as 



ip = C123 



dX $,dv 
dt dr 



C123 



dA + + ^ _ _2_ 

dr dr I r 2 



(83) 



Equations (f8"2"|) and (j53"|) can then be combined to elimi- 
nate any explicit v terms: 



8X 
dt 



r 
a 



-C123 



£123 (&± (OA 

2 \dr 2 \dr 
dAdt> | /c 123 

V C14 



dr 



dip 
dr 



C123 - c 2 + 1 



IdA 

r dr 



c 2 



1 d<S> 

r <9r 



= (84) 



where we have defined the constant a to be 
a=(ci+c3 + l)(ci 23 + 2c 2 -2). 



(85) 



We can then use this equation and (|81|) to write down a 
single second-order time-evolution equation for tp: 



e 2A-2$ r 2_£i£ 



d 2 ip 



r 
a 



-C123 



+ C14 at 2 
Cl23 ^ + ^ 

9A9£ + /ci23 _ 
9r dr \ C14 

/ C123 
V C14 



9$ 



ci + c 3 + 1 



C 2 



1 <9$ 

r <9r 



dtp_" 
dr 

IdA 

r dr 

= 0. (86) 



Moreover, using (|T8[) to eliminate the d(p/dr term from 
([70]) . we find that 



^-2ci 23 e- 2 *(A 2 V'i-Ai^ 2 ) (87) 
and the reduced symplectic form can be written as 

C123C14 



-47T- 



C14 Jo 



drr 4 e 3A ~ 3 * 



cVi , dtp 2 , 



The form W a p thus defined may be positive or neg- 
ative definite, depending on the signs of ci 23 and ci 4 
(recall that we are assuming that Ci 23 and C14 are non- 
vanishing); however, it is never indefinite. Thus, the sym- 
plectic form is in fact of the required form (|16p . and the 



11 



equations of motion can be put in the form (fTJ). We can and the numerator is 
therefore write down our variational principle of the form 
([3|); the denominator is 

(V>» = -4^^Ii/ drr 4 e 3A - 3 V (89) 

2 + c 14 Jo 

I 



C123 



(</>, T if;) = -4tt— 



j 4 A-<t> 

dr re 



C123 



<9i/> 
9r 



C123 



OAd$ _ /C123 _ 



<9r dr 



C2 



1 <9$ 

r <9r 



C123 
C14 



Ci + C 3 



ISA 
r <9r 



(90) 



C. Discussion 

We can now examine the properties of this variational 
principle to determine the stability of Einstein-aether the- 
ory. The simplest case to examine is that of flat space. 
Let us denote the coefficient of ip 2 in (j9"0|) by Z(r), i.e., 



Z(r) = ci23 



dr dr 



<9A<9£ _ / ci23 

\ C14 
/ C123 
V Cl4 



C2 



r or 



+ Cl + c 3 + 1 



r or 



In the case of fiat spacetime, Z(r) vanishes, and we are 
left with 



C123 



(2 + c 14 )/o 00 d- 4 (f) i 



J Q dr r 4 i/> 2 



(92) 



Since both integrands are strictly positive, we conclude 
that flat space is stable to spherically symmetric pertur- 
bations in Einstein-aether theory if and only if 



Cl23(2 + C 14 ) 



ci 4 (ci + c 3 + l)(ci2 3 + 2c 2 - 2) 



> 0. 



(93) 



This is the same condition found in 0] for the stability of 
the "trace" wave mode (modulo the aforementioned sign 
differences) . The other stability conditions found in Q - 
those corresponding to the "transverse traceless metric" 
and "transverse aether" wave modes — are incompatible 
with spherical symmetry, and thus our analysis cannot 
reproduce these conditions. 

More broadly, we can also analyze the stability of the 
general spherically symmetric solutions described in [24| . 
These (exact) solutions are described in terms of a pa- 
rameter Y: 



r(Y) 



Y-Y- (Y -Y 



Y 



Y -Y+ 



(94) 



$(Y) 



Y + 



2(2 + Y+ 



■In 



i - 

i - Y/Y+ 



A(y) = i in (-^(y-y + )(y-y_; 

where the constants Y± are given by 



(95) 



(96) 



(97) 



and r m i n is a constant of integration related to the mass 
M via 



2M 

~-y1 



(_\ _ y+ )(l+y+)/(2+Y+) 



We can then obtain Z(Y) by writing out (|9~lj) in terms of 
these functions of Y (noting that, for example, dA/dr = 
(dA/dY)/(dr/dY)), and then plot Z and r parametri- 
cally. The resulting function is shown in Figure [1] for 
C123 = ±ci4 and C14 = —0.1. 

The first thing we see is that in the asymptotic region, 
the sign of Z(r) is determined by the sign of C123. To 
quantify this, we can obtain an asymptotic expansion of 
Z(r) as r — * 00, noting that r — » 00 as Y — > 0. Doing 
this, we find that to leading order in M/r, 



Z{r) 



M 

Cl23^3" + 



(99) 



Thus, the coefficient of the ip 2 term in (|90[) is, to leading 
order in the asymptotic region, the same sign as that 
of the (dip/dr) 2 term. We can therefore conclude that 
for a test function in the "asymptotic region" of one of 
these solutions (i.e., r ^> M), the conditions on the c;'s 
are the same as those in flat spacetime. 8 For a realistic 



Note that if the leading-order coefficient in 11991 had been differ- 
ent from that of the (^) 2 term in (1906 , we would have obtained 
a new constraint on the Ci coefficients; the matching of these two 
coefficients seems to be coincidental. 
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FIG. 1: Representative plots of Z(r) versus r for the solu- 
tions described in [24|]. The solid and dashed lines corre- 
spond to C123 = ±0.1, respectively; in both cases, C2 = and 
C14 = —0.1. Note that the choice of parameters leading to 
the dashed plot would lead to instability in the asymptotic 
region. 



spacetime, of course, there will be some ball of matter 
in the central region, and the region where the vacuum 
solution holds will be precisely the region where r 3> M; 
thus, we can conclude that for a normal star in Einstein- 
aether theory, the exterior is stable if and only if 
holds. (Henceforth, we will assume that the Cj coefficients 
have been chosen with this constraint in mind, unless 
otherwise specified.) 

We also note that as r approaches r m - m , Z(r) diverges 
negatively. One might ask whether a test function in this 
region (in a spherically symmetric spacetime surrounding 
a compact object, say) could lead to an instability. We 
investigated this question numerically using our varia- 
tional principle; however, our results for (ip, Tip) were 
positive for all test functions ip that we tried. Roughly 
speaking, the derivative terms in (|89[) always won out 
over the effects of the negative Z(r). 9 This is, of course, 
far from a definitive proof of the positivity of ([50)1 ; and 
it should be emphasized that the above analysis has not 
considered the effects of matter on the stability of such 
solutions. Nevertheless, the above results are at least in- 
dicative that the spherically symmetric vacuum solutions 
Einstein-aether theory do not possess any serious stability 
problems. 



9 We have been assuming a "static aether" throughout this work; 
however, there also exist solutions in which the aether vector 
u a is not globally aligned with the time-translation vector field 
[25l |. Thus, even supposing that such "extremely compact stars" 
are phenomenologically realistic, and that the exterior solutions 
for such objects are unstable, these "non-static aether" solutions 
might still be physically viable. 



D. The case of C123 = 

The above analysis assumed that C123 was non- 
vanishing; however, the action originally considered by 
Jacobson and Mattingly [§[ used a "Maxwellian" action, 
i.e., the kinetic terms for the aether field were of the form 
F a bF ab , corresponding to c 3 = — C1 and c 2 = 0. This case 
is therefore of some interest. 

In the preceding analysis, the perturbational equa- 
tions of motion (|79|) and (|80|) were obtained without any 
assumptions concerning C123, as was the preconstraint 
equation (|75|) . We can therefore simply set C123 to zero 
in these equations and use ([75)) to solve for d<p/dr, as we 
did in the general case. After application of the back- 
ground equations of motion, there result the equations 



2 



1 A + e 



2A-* ' 



■(ca 



„ dv 



and 



-(l-ca)— -Z v = 0, 
r at 



(100) 



(101) 



where we have defined Zq in terms of the background 
fields: 



Z = — (—) - ~C 2 (— + — 

2 V dr J r \ dr dr 



2dA 

r dr 



1 



,2A 



1 



(102) 



We can then combine these two equations to make a sin- 
gle second-order equation for v. 



p 2A-2<i> gi4(c2 ~ I) 2 d 2 v 



C14 



<9t 2 



(103) 



Since there are no radial derivatives to contend with, the 
solutions to this equation are all of the form 



v(r,t) = /(r)exp 



±\ Mr) 



Cl4 



,*-A 



c 2 - lit 



(104) 

where f(r) is an arbitrary function of r. 

Stability of these solutions will thus depend on the 
sign of the quantity + l)Z (r). We can plot Z (r) 
parametrically, as was done in the general case; the re- 
sulting function, shown in Figure [2l falls off rapidly in 
the asymptotic region. We can also perform an asymp- 
totic expansion similar to that done in the general case 
to find the large- r behaviour of Z (r); the result is that 



M 2 

Z (r) w (1 - c 2 )ci4^- 



(105) 



We conclude that in the C123 = case, the spherically 
symmetric static solutions of Einstein-aether theory are 
unstable unless 



(c 2 -l)(ci4 + 2) >0, 



(106) 
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FIG. 2: Plot of (^2_ + \)Z (r), relevant for Einstein-aether 
theory with C123 = d. In this plot, we have chosen cu — —0.1 
and C2 = 0.1. 



i.e., unless either C2 > 1 and C14 > —2, or C2 < 1 
and C14 < —2. Note that this excludes the limit of 
"small Cj" often considered in works such as Q- In the 
"Maxwellian" case, this instability is likely related to the 
non-boundedness of the Hamiltonian 2jj; however, our 
analysis also suggests that non-standard kinetic terms 
can in fact stabilize the theory (at least in the time- 
like aether case), contrary to the arguments made in that 
work. 

Finally, we note that even if (| 1 06[) holds, the perturba- 
tional solutions of the form (|104p may still be problem- 
atic, as the radial gradients of v will grow linearly with 
time. To quantify what magnitude of gradients would be 
acceptable, we note that in the C123 = case, there is 
only one dynamical degree of freedom; thus, v uniquely 
determines A and d<j>/dr. Setting C123 = in ([78]) and 
(|75|) allows us to solve these equations for A and d(j>/dr; 
in particular, 



dr 



J2A 



-cue 



2A 



2 3$ 

C14 — 

r or 

2~ 



C14 



dr 



2 

-c 2 



„2A-* 



C14 



d® 

dr 



dv 

at' 



which simplifies in a nearly-flat spacetime to be 

C14 + 2c2 dv 



dr 



Ci 4 + 2 dt' 



(107) 



(108) 



Applying this to our solution (|103[) . we find that for suf- 
ficiently large t 



d 2 



dr 2 



2y/(c 2 -l) 3 (ci4 + 2c 2 ) 



M 



dtp 



dr 



t 



(109) 



We can estimate typical scales for d<j>/dr by looking at 
other sources, such as planets. The scale of perturbations 
due to the Earth's gravitational field (at its surface) is 
given by d<j>/dr « M^/r^, where M e and r® are the 
mass and radius of the Earth, respectively. Thus, the 



perturbation to d 2 cj)/dr 2 near the Earth's surface will be 
of the order 



d 2 



dr 2 



2V( C2 -l) 3 ( C i4 + 2c 2 )^ 

It Tm 



t 



(110) 



where R is the radius of Earth's orbit and M Q is the mass 
of the Sun. This enhancement to d 2 4>/dr 2 will lead to an 
observable change in the tidal effects due to the Sun's 
gravity; demanding that these remain small relative to 
the normal Newtonian tidal effects, i.e., d 2 cj)/dr 2 <c 
d 2 <$>/dr 2 2M Q /i? 3 , we then have 



v/(c2-l) 3 (ci4 + 2c 2 ) « 



1.9 x lo- 



an) 



for t w 5 x 10 9 years (the approximate age of the Earth.) 
We see, therefore, that the requirement that perturba- 
tional tidal effects remain small severely constrains our 
choices of C14 and C2. 



V. TeVeS 
A. Theory 

TeVeS (short for "Tensor- Vector-Scalar") is a modi- 
fied gravity theory proposed by Bekenstein [8j in an at- 
tempt to create a fully covariant theory of Milgrom's 
Modified Newtonian Dynamics (MOND). The fields 
present in this theory consist of the metric; a vector field 
u a which (as in Einstein-aether theory) is constrained by 
a Lagrange multiplier Q to be unit and timelike; and two 
scalar fields, a and a. The Lagrangian four-form is 



where C g is the usual Einstein-Hilbcrt action, 



C s is the "scalar part" of the action, 



(112) 



(113) 



C s = -^a 2 (g ab -u a u b )V a aV b a~r 2 a 4 F(ka 2 ), (114) 

with k and £ positive constants of the theory (with 
"length dimensions" zero and one, respectively), and 
F(x) a free function; C v is the vector part of the action, 

Cv = -^-F ab F ab + Q(u a u a + 1), (115) 

dill 

with K a positive dimcnsionless constant and F ab — 
V a iifc — Vf,u a ; and C m the matter action, non- minimally 
coupled to the metric: 

C m e - C mat [A,e 2a g ab + 2u a u b sinh(2a)] (116) 

where jC mat [^4, g ab ] would be the minimally coupled mat- 
ter Lagrangian for the matter fields A. Note that if we 
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ignore the scalar and matter portions of the Lagrangian, 
this Lagrangian is the same as that for Einstein-aether 
theory (|64|). with c\ = — C3 = — K/16n and c 2 = C4 = 0. 
In the present work, we will work exclusively with the 
"vacuum" (£ m = 0) theory. 

Taking the variation of (1112j) to obtain the equations 
of motion and the symplectic current, we find that 



-*-A 



and 



£a = 



dr 



a 2 r 2 e*- A 



da 
dr 



(120c) 



(120d) 



5C = ((£ G )abSg ab + (£ u )aSu a + £ a 5a 



-£Ja + £ Q 5Q + V a 8 a )e, (117) 



where 



{£G)ab = 77— G ab - \a 2 \7 a aS7 b a - Qu a u b 

+ A- ^,„V7 i)c - F ac F b c + ±g ab F cd F cd ^J 

+ \g ab Qa 2 (V c aV c a - a 2 ) + V V 4 F(^ 2 )) , 

(118a) 



where we have defined y(x) 
0) and used the equation for Q 



xF(x) — \x 2 F'{x) (as in 



Q 



K 

loTf' 



-2A 



<9 2 $ d§ fdh 2 
9r 2 dr \ dr r 



(121) 



to simplify. (As in Einstein-aether theory, the background 
equation (£ u ) a = is satisfied trivially by a static aether.) 

The symplectic form for the theory can now be calcu- 
lated from (|119p . The result will be essentially the same 
as that in Einstein- aether theory (with the appropriate 
values of the Cj's), with added terms stemming from the 
variations of C s : 



{£ u ) a = a 2 aV a a + ^-V b F ba + 2Qu a , (118b) 



£ a = V a (a 2 (V a a-u a a)), 



7 Ein 



(122) 



where Wg in is given by ([55)1. lj" cc is given by (|B1[) with 



(118c) c i = — c 3 = — K/IQtt and C2 = C4 = 0, and 



V"aV„a - a 2 



1 



+r 2 a 2 [F(ka 2 ) + -ka 2 F\ka 2 
and 

Q a = Okn ~ CT 2 (V a a - u a a)5a 

+ §- (F b a Su b + F\u c Sg bc ) 



(118d) 



(119) 



In the above, we have defined a = M a V a a. The remain- 
ing equation £q = is identical to that in Einstein-aether 
theory, (|67c|) . For a static solution, in our usual gauge, 
the background equations of motion become 



(£ 



Gltt 



„2$-2A 



1D7T \r or r z 




<9 2 $ 5$ fdk 2 
dr 2 dr \ dr r 




, (120a) 



(V"o - 2^-i 9bc ^ c 



+ <5i 5 afc V 6 a + ( 3 Qb - u a u b )Vf,(5ia 



2u {a 5 1 u b) V b a 



5 2 a. (123) 



B. Applying the variational formalism 

The next step is to write out the ^-component of uj a 
in terms of the perturbational variables. We will take 
our metric perturbations to have the usual form, and our 
vector perturbation to be of the same form as was used 
for Einstein-aether theory (fTTj) . For the two scalar fields, 
we define 5a = (3 and 5a = t. Calculating w' in terms 
of these perturbational variables, and using the results of 
([75]). we find that 



at 



1 * 

+e v\ 



k (a® 

^Vdr~ X2 



dr 



(124) 



re ^ 1/2 9$ 1 2A A if /9$ 

(£ G )rr = T^ To a( e ~ *) + 



167r \r dr 

1 2 / da 
~4 a (dr 



16tt \ d 



e2 V' j4F(fcff2) ' (12 ° b) 



We now turn to the question of the constraints. The 
tensor C ab is again given by C ab = {£c)ab + u a (£ u ) b , 
and thus to first order we have 5C rt = {£c)rt + 3u r (£ u ) t 
(since u r = in the background.) Calculating 5C r t in 
terms of our perturbational variables, we find that the 
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preconstraint equation is 



F = r 2 e*" A 



, da 



1 



Or 



.9$ 



CT -5-P-7T ~+K^- A 



K 

"ST 



8tt \r 
,dv 
dt 



,2A-* 



Or 
K d(j> 
8ir dr 



0. (125) 



In principle, we could now use this equation, together 
with the equation 



„2A / a cVr„_2\ 1 



+ .(_L^ W) _(*' 



A 



„2A-* 



2 9a 9/3 1 /2 9$ 
9r dr 8tt \ r dr I dr 



K 9$ dv 
8tt 9r ~dt 

(126) 



to obtain equations for A and d<p/dr in terms of the "mat- 
ter" variables v, (3, and r. However, it is simpler to pur- 
sue a similar tactic to the one we used in the reduction 
of Einstein- aether theory. To wit, the r-component of 
(6£ u ) a is 



(S£ u ) r = \^>e 2A Q + <7 2 
K /9$9A 




-3> 2 

v + e a 



2A-$ 



d 2 v d 2 



dadf3_ 
dr dt 

= 0. (127) 



iir \ dr dt dt 2 drdt J 

We can use (|125p to simplify this equation; the result is 

2\ 



^(J_ X _ 2 A ] (128) 

dt \ Airr dr ' ' y ' 



We then define the new variable x a s 

X= 7-A-2a 2 — /3. 
47rr or 



(129) 



The evolution equation for f3, meanwhile, is given by 
the equation 5£ a = 0; in terms of the perturbational 
fields, it is 



2e z 



d 2 (3 


d 2 /3 


(Ida 


dt 2 + 


dr 2 


\a dr 


da 


fd4> 


9A\ 


dr 


[dr 


dr J 



1 da da 
a 2 dr dr 



„2A-$ 



9$ 9A 
dr dr 
2 da dr 
a dr dr 
da dv 



2\ 9/3 



dr dt 



r J dr 



0. (130) 



Finally, the field r, being the perturbation of the auxil- 
iary field tr, can be solved for algebraically in the equation 
8£ a = 0: 



r 2 ay'(ka 2 )T + e 




da 9/3 
dr dr 



0. (131) 



We can thus follow the following procedure to write the 
evolution equations solely in terms of (3 and x : first, 
we use the preconstraint equation (|125| to eliminate the 
combination d(f>/dr + e 2A ~*9w/9i in favour of A and (3; 
next, we use (| 13 1[) to eliminate r and its spatial deriva- 
tives; and finally, we use the definition of x (|129|) to elim- 
inate A and its derivatives in favour of (3 and x- Applying 
our procedure to (|126p results in an equation depending 
on dv/dt as well as (3, d(3/dr, and xi combining this 
equation with (|128p then yields a second-order evolution 
equation for x of the form 



d 2 X w , v 9/3 



V 4 /3. 



(132) 



where the coefficients V$ and W are dependent on the 
background fields. Similarly, applying this procedure to 
()130|) . and using (|126|) to eliminate a resulting dv/dt 
term, we obtain a second-order equation for (3: 

d 2 (3 d 2 9/3 9 X 

■^2 = W i^ 2 -+ W 2^:+^ + V 1 — + V 2X , (133) 

where, again, the Wi's and Vj's depend on the background 
fields. Note that (1132|) has no dependence on the spatial 
derivatives of x, only upon \ itself. 

We also need to reduce the symplectic form and ex- 
press it in terms of (3 and x- Applying the preconstraint 
equation (| 125[) to (|124p . and using the definition of \ 
(fT29]) along with (fl28l) . we find that 



fl = 4vr / drr 2 e A "* 



1 ^Xl „ 2 9 ft«' 



where 



[1h2], (134) 



(135) 



We can also apply the background equations of motion, 
along with the Equation (|121[) , to this quantity to obtain 



H 



1 /9A 9$ 



47rr \ dr dr 



— + — - 2a' — 



dc 



dr 



(136) 



For a variational principle to exist for this theory, the 
form W a/ 3 defined by Q134|) must be positive definite. 
In our case, this means that the coefficients of both 
(dxi/dt) X 2 and (9/3i/9t)/3 2 in the integrand of (fTM)) 
must be positive for the background solution about which 
we are perturbing. While the coefficient for the latter 
term is obviously always positive, the situation for the 
former coefficient (namely, H) is not so clear. To address 
this issue, we need to know the properties of the spher- 
ically symmetric static background solutions of TeVeS. 
These solutions (with a "static aether") are described in 
[28| . In our gauge, they are most simply described in 
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H(r) 
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FIG. 3: Plot of H(r) versus r for TeVeS. In this plot, k and 
K are both chosen to be 10 -2 , and m s = z g = 1. 



terms of a parameter z 



r(z) = 



10 



z + z c 



-Z g /4:Z C 



Zg/iZo 



(z 2 + z c 2 ) 



a(z) 



km s 
8nz c 



In 



(137) 



(138) 



(139) 



(140) 




where z c , z g , m s , and a c are constants of integration. 
The first three of these are related by 



(141) 



while a c "sets the value of a at oo." The constant m s , 
which can be thought of as the "scalar charge" of the 
star, is defined by an integral over the central mass dis- 
tribution Q; for a perfect fluid with p + 3P > 0, m s is 
non-negative. Finally, the constant z g is also defined in 
terms of an integral over the central matter distribution 

1- 

Plotting H(r) parametrically (Figure [3]), we see that 
this function is strictly negative. In fact, in the z — > oo 
limit, we have r/z = 1 + 0(z c /z), and so we can take 
r w z to a good approximation. Calculating H(r) in 
terms of $(r) and A(r), we find that as r — > oo, 



H{r) 



1 

'16tt 



k ~ml + f zl 



1 



(142) 



which is negative for any positive choice of k and K. 
Therefore, it is not possible to straightforwardly apply 
the variational principle to TeVeS, since the quadratic 
form defined by W Q ^ is indefinite and thus cannot used 
as an inner product. 11 



C. WKB analysis 

While we cannot derive a variational principle with 
TeVeS, the fact that we have been able to "reduce" the 
equations of motion to an unconstrained form still allows 
us to analyse the stability of its spherically symmetric 
solutions. Let us consider a WKB ansatz, of the form 



(3{r) 
x(r) 



i(u){r)t+nr) 



'Mr) 
Mr) 



(143) 



with k very large compared to the scale of variation of 
the background functions Ui, Vi, and W. We will further 
choose the the functions fi(r) and oj(r) are chosen to be 
"slowly varying" relative to the scale defined by k, i.e. 



1 dft 



< At, 



1 duj 



< K. 



fi(r) dr ^(r) dr 

Under this assumption, we will then have 
d 



(144) 



fe*W Hw ' /i(r)) w iK,e l ^ ( - r)t+Kr) j l {r). (145) 



Now let us apply the time-evolution operator T im- 
plicitly defined by (fT3"3"|) and (fl^|) to our ansatz (fT¥5|) . 
We see that for sufficiently large k, the highest-derivative 
terms will dominate the lower-derivative terms. Thus, to 
a good approximation we will have 



T 



7?" 




X. 





—K 2 Ui(r) ZKVi(r) 
inV 3 {r) W(r) 



'Mr) 
Mr) 



J(u(r)t+K,r) 



(146) 

Then, in the limit of large k, our ansatz (|143|) will be an 
approximate eigenvector of T if there exist a fp(r) and 
f x (r) such that 



-u 2 (r) 



'Mr) 




Mr)_ 





-K 2 U\{r) inVi(r) 
kV 3 (r) W{r) 



Mr) 
Mr) 



(147) 



In other words, in the limit of large K, the problem of 
finding modes of T is a simple two-dimensional eigenvalue 
problem where the eigenvalues are functions of r. In this 
limit, the eigenvalues of this matrix are (to leading order 
in k) 



w 2 (r) w <^ K 2 Wi(r),-W(r) + 



Vi(r)V 3 {r) 
Ui(r) 



(148) 



This z is the "radial coordinate" in isotropic spherical coordi- 
nates, as used in |2S| . 



Note that this also implies that the "perturbational Hamilto- 
nian" of TeVeS, as defined in equation (54) of has an indefi- 
nite kinetic term. 
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It is not difficult to verify that for sufficiently large k, our 
assumptions for the ansatz (|144[) are satisfied. 

It remains to write out the functions Ui(r), Vi(r), 
V3(r), and W(r) in terms of the background functions. 
These can be shown to be: 

Wi = - e 2 *- 2k I 1 + e -2A 2£2 (If) 2 ) / 1491 

v, - -i>-^ (l + e -^^ML\ , m) 

V 3 = 16^re 2 *- 2A a 2 ^ 
or 

2 /3a \ 1 fdA M . 

and 

W = ,V— (|-!- 8 ^(^) 2 ) 

x { 8 ™ uo -n^ + ^)J- (152) 

This implies that the eigenvalues of the matrix in (I147|) 
are 

k- 2 / 2£ 2 f^ 2 \ 

" ! W-T«"-"(i + «- 2A ^Pyj ("») 

and 

- 2 W « e 2 *- 2A (| - l) 

/ 1 /<9A <9$\ 2 fda\ 2 \ , . 

We can see that this first eigenvalue (|153|) is always 
positive as long as y'(x) > 0; indeed, the choice of y(x) 
made in [1[ does satisfy this inequality. The second eigen- 
value (|154p . however, is just 

u 2 (r) = 4^e 2 *- 2A (J^ - lj H{r). (155) 

Since H(r) is always negative, we conclude that this sec- 
ond mode is unstable for < K < 2. Further, for a 
spherically symmetric solution outside a Newtonian star, 
the approximation (|142p is valid; thus, to lowest non- 
vanishing order in r _1 , we have 

- 2 «--j(|-l)(^ + f, fl 2 )l (156) 



D. Discussion 

The above result for Lo 2 (r) implies that for stability of 
Newtonian solutions, we must have K > 2, ruling out 
the range of parameters originally considered in |8J] . To 
estimate the time scale of the instability when K < 2, 
we first note that unless k and K are much larger than 
the ratio of the star's radius R to its Schwarzchild radius 
m g , the parameters m g , m s , and z g are all approximately 
equal, differing by terms of 0(m g /R) Q. In particular, 
if k and K are of the same order and K 2, we find that 
the timescale of this instability will be on the order of 
10 6 seconds — approximately two weeks — for points near 
the surface of the Sun. 

Note also that the effective gravitational potential for 
nonrelativistic motion is 

$ = E$ N + a, (157) 

where 

Z = e- 2a -(l-K/2)-\ (158) 

is the Newtonian gravitational potential, and a c is 
the asymptotic value of a, determined by the cosmologi- 
cal boundary conditions [29] . Expanding the asymptotic 
behaviour of <&jv and a about infinity, we find that 

MfTi + 5) ? (15<J) 

where in the weak-field limit, we have 




(160) 



and m s » e~ 2ac m g . The coefficient in parentheses in 
(|159p can be thought of as renormalizing the gravita- 
tional constant G by some overall factor; the first term 
in parentheses comes from the usual perturbations to the 
"physical metric" component gtt, while the second term 
comes from variation of the scalar field a. However, if 
K > 2, this would make the metric contribution to the 
"effective gravitational constant" negative. This could 
be remedied by choosing k > 8ir/(K — 2). However, it 
is unclear whether the theory, in this parameter regime, 
would still be experimentally viable; for example, Big 
Bang nucleosynthesis constraints require that k < 0.75, 
and the CMB power spectrum also constrains the viable 
regions of parameter space [30]. 12 



We also note that if the gravitational effects of the scalar are 
dominant over those of the tensor, the solutions of TeVeS would 
greatly resemble those of the so-called "stratified theories" with 
time-orthogonal space slices [3l| . As these are ruled out exper- 
imentally via geophysical experiments, we could not make k or 
K too large without running afoul of these experimental con- 
straints. 
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We also note that the eigenvector corresponding to the 
unstable mode of T will satisfy (to leading order) 



Vi 
Ui 



fx 







(161) 



or, in the limit of large n, fp k, 0. In other words, this un- 
stable mode should manifest itself in growth of the radial 
component of the vector field u a rather than growth of 
the scalar a (cf. (|128ll 1. An instability of the vector field 
was also found in [32] by Dodelson and Liguori. How- 
ever, it seems unlikely that this is the same instability 
for three reasons. First, the instability found in was 
found in a cosmological context, not a Newtonian-gravity 
context; since the cosmological solutions of TeVeS are in 
a very real sense separate from the Newtonian solutions 



(existing on two different branches of the function y(x)), 
it is difficult to draw a direct correspondence between the 
stability properties of these two types of solutions. Sec- 
ond, the instability found in [32] manifests itself only in 
the limit of a matter-dominated Universe; the instabil- 
ity we have found exists in vacuo. Third, Dodelson and 
Liguori's instability requires K to be sufficiently small 
relative to k, while our instability is present for all k > 
and < K < 2. 

Finally, it is perhaps notable that the vector field u a in 
TeVeS has "Maxwellian" kinetic terms, which we found 
in Section IIVI to be unstable in the context of Einstein- 
aether theory. It is possible, then, that this instability 
could be cured via a more general kinetic term for the 
vector field. 
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APPENDIX A: COEFFICIENTS FOR THE f(R) VARIATIONAL PRINCIPLE 



The quantity (ip,Tip) for f(R) gravity can be put in the form 



(ip, T tp) = 4?r / dr 



Ci 



where the C, coefficients are given by 



db 



d = 6e*- A S- 2 /'(a), 
C 2 =87rr 2 e A +VV 2 , 
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(A6) 
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and 



C 6 = Stt^+VS- 3 { S 2 



{g') 2 d$ (d§ 2 



g" dr \dr r 



Q v [-2 + ~ 



dA /<9$ 2 



r dr V dr J dr \ dr r 



3 // 2 



e 2A ^V -(/(a)-16vr e ) + 8V- , 
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APPENDIX B: SYMPLECTIC CURRENT FOR EINSTEIN- iETHER THEORY 

The symplectic current for vacuum Einstein- aether theory is given by uo a — Wg in + wj ec , where Wg in is the standard 
Einstein symplectic form (given by (|36[) ). and 

< c = (M (1) ) a bc d ef 6 2 g bc V d 6 ig ef + (M {2) ) a bcde S 2 g bc S ig de + (M (3) )V e S 2 g bc V ' d 6 lU e + (M (4) ) a bcd S 2 g bc S lU d 

+ (M {5) ) a b c de 5 2 u b V c 6 ig de + (M (6) ) a b c d 6 2 u b V c S lU d + (M (7) ) a bc 6 2 u b S lU c - [1 «-> 2], (Bl) 

where tensors are given by 

/ 7i ,t \a d ( ad ca cd , a d \ , ^ ad 

(-M(i)j 6c e/ = ci I g g be u c u f - d e d b u c u f + -u u g be g c f J + ^c 2 u u g bc g e f 

+ c 3 (s a e S d b u c u f - g ad g be u c u f + ^u a u d g be g cf ^j - c 4 ^2u a u d g be - S a e u d u b - u a u e 5 d b + ]^g ad u b u^j u c u f , (B2) 



(M {2) ) a bcde = (c 3 - a) (v b u a g ce u d + 5 a d u c V e u b + u a g bd V c u e + ^g de u c (\7 b u a - W a u b )^j 



+ 2 ( c i + c 3) ua V b u c g de - c 4 (u a ii c - -ii a u c ) u b y, : , . ( Bo ) 



(-^(3)) a 6c d e = ci (if^^iiic - g ad g be u c - u a S d b g ce ) - c 2 u a g bc 5 d e 

+ C3 (g ad 9beU c - S a e S d b u c - u a S d b g ce ) - c 4 {S a e u c - 2u a g ce ) u b u d , (B4) 

(M( 4) ) a fccd = ci {W b u a g cd + W a u b g cd + V a u d g bc - 8 a d W b u c - 2<5 a c V fc u ci ) 

+ c 3 (V a u b 3 cd - V b u a g cd - S a d V c u b + V d u a g bc ) 

- c 4 (2u a u b g cd + u b u c V d u a ~ 25 a d u b u c - 2u a u b V d u c + u a g bc u d ) , (B5) 



(M (5) )V de - K a db c u e - K ac bd u e - K a dbe u c , 



(M (6) ) a b c d = 2K ac bd , 



(B6) 
(B7) 



and 



(M (7) )\ c = -2c 4 {8 a c u b + u a V c u b ). 



(B8) 
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